home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / REDUC.z / REDUC
Text File  |  1996-03-14  |  3KB  |  67 lines
  1.  
  2.  
  3.  
  4. ____RRRREEEEDDDDUUUUCCCC((((3333FFFF))))                                                          ____RRRREEEEDDDDUUUUCCCC((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      REDUC, SREDUC   -  EISPACK routine.  This subroutine reduces the
  10.      generalized SYMMETRIC eigenproblem Ax=(LAMBDA)Bx, where B is POSITIVE
  11.      DEFINITE, to the standard symmetric eigenproblem using the Cholesky
  12.      factorization of B.
  13.  
  14.  
  15. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  16.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  rrrreeeedddduuuucccc((((nnnnmmmm,,,, nnnn,,,, aaaa,,,, bbbb,,,, ddddllll,,,, iiiieeeerrrrrrrr))))
  17.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
  18.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn aaaa((((nnnnmmmm,,,,nnnn)))),,,, bbbb((((nnnnmmmm,,,,nnnn)))),,,, ddddllll((((nnnn))))
  19.  
  20.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssrrrreeeedddduuuucccc((((nnnnmmmm,,,, nnnn,,,, aaaa,,,, bbbb,,,, ddddllll,,,, iiiieeeerrrrrrrr))))
  21.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, iiiieeeerrrrrrrr
  22.           rrrreeeeaaaallll             aaaa((((nnnnmmmm,,,,nnnn)))),,,, bbbb((((nnnnmmmm,,,,nnnn)))),,,, ddddllll((((nnnn))))
  23.  
  24.  
  25.  
  26. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  27.      On Input
  28.  
  29.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  30.      as declared in the calling program dimension statement.
  31.  
  32.      NNNN is the order of the matrices A and B.  If the Cholesky factor L of B is
  33.      already available, N should be prefixed with a minus sign.
  34.  
  35.      AAAA and B contain the real symmetric input matrices.  Only the full upper
  36.      triangles of the matrices need be supplied.  If N is negative, the strict
  37.      lower triangle of B contains, instead, the strict lower triangle of its
  38.      Cholesky factor L.
  39.  
  40.      DDDDLLLL contains, if N is negative, the diagonal elements of L. On Output
  41.  
  42.      AAAA contains in its full lower triangle the full lower triangle of the
  43.      symmetric matrix derived from the reduction to the standard form.  The
  44.      strict upper triangle of A is unaltered.
  45.  
  46.      BBBB contains in its strict lower triangle the strict lower triangle of its
  47.      Cholesky factor L.  The full upper triangle of B is unaltered.
  48.  
  49.      DDDDLLLL contains the diagonal elements of L.
  50.  
  51.      IIIIEEEERRRRRRRR is set to Zero       for normal return, 7*N+1      if B is not
  52.      positive definite.  Questions and comments should be directed to B. S.
  53.      Garbow, APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.